Universal free energy correction for the two-dimensional one-component plasma 
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r-| ! I. INTRODUCTION 

o 

Two-dimensional Coulomb systems are models which have attracted some attention. On a two-dimensional mani- 
fold, a Coulomb system is a system made of particles interacting through the corresponding Coulomb potential plus 
perhaps some short-range interaction. In a plane, the Coulomb interaction energy of two particles of charges q and 
^ ■ q' , separated by a distance r, is defined as —qq'\n{r/L), where L is some (irrelevant) length. 

c/2 [ Some time ago, it has been shown that the free energy of such systems has a universal finite-size correctionjQ very 
similar (except for its sign) to the one which occurs in a system with short-range forces at a critical point 0: for a 
finite Coulomb system of characteristic size R, the free energy F has the large- i? behaviour 
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O ] where /3 is the inverse temperature. A and B are non-universal constants describing the bulk and boundary contri- 
butions, respectively. (x/6)lni? is the universal correction, depending only on the Euler number x which describes 
the topology of the manifold on which the system lives. However, the general derivation of (|l]) had some heuristic 
features. The purpose of the present paper is to check (|l|) in the special case of a one-component plasma on the 
K*" ■ surface of a sphere, by a different method. 

The one-component plasma is a system made of one species of point-particles of charge g in a uniform neu- 
■ tralizing background. On a sphere of radius R, the interaction between two particles can be chosen |^,^ as 
' — ln[(2_R/L) sin(?/'/2)], where ip is the angular distance (seen from the sphere centre) between the two particles. 
There are also particle-background and background-background interactions. A dimensionless coupling constant is 

O r = /3<72. 

, A sphere has no boundaries and its Euler number is x = 2. Furthermore, for a given particle density, R^ is 
' proportional to the number of particles N . Thus, expansion (^ becomes 

^ pF = CN+lliiN + ---, (2) 

^ , where C is a constant. The model is exactly solvable Q when F = 2 and it can be checked that (||) is obeyed in 
O ■ that case. Also, exact calculations Q for finite values of iV at F = 4 and F = 6 are well fitted by (||). 

' The present derivation of (||) relies on a recent result ^ about the direct correlation function c(r) of the plane 
one-component plasma. By a diagrammatic analysis, it has been shown in ref. |^ that the second moment of the 
short-range part Cg^(r) has the simple value 
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A remarkable feature of is its universality, in the sense that it is independent of the coupling constant T. It will 
now be shown how (H) leads to (||) . More specifically, we show how (^) implies the finite-size correction to the chemical 
potential fi = dF/dN: 

= /3a/oc + ^ + . . . (4) 



This derivation bears some similarity with another one about the two-component plasma [TTUg 



II. DENSITY FUNCTIONAL THEORY APPROACH 



We consider the OCP of average density Us on the sphere of radius R (with a corresponding number of particles 
N — A-kR^ Us). Introducing the stereographic projection of the sphere onto the plane V tangent to its south pole (see 
figure 1), we map the homogeneous OCP on the sphere onto a modified inhomogeneous plasma on the plane, with 
local particle density 



n(r) 



r 

IR2 / 



V 



(5) 



In terms of planar coordinates ri and r2, the interaction potential between two particles on the sphere with angular 
distance V'12 can be written as the sum of the planar two dimensional Coulomb potential Vp{r, r') = —q^ ln[|r — r'\/L] 
and one-body terms since: 



-In 



2R 



sm 



2 



= - In 



|ri -r2| 



L 



ilnd 



A.) 



ilnd 



(6) 



The two one-body terms appearing on the right hand side of eq. as well as the metric, create a central potential 
and we can consider the projected planar system as a OCP interacting through the standard pair potential Vp, in 
an external one-body central potential (r). The latter acts as a confining mechanism ensuring the proper density 
given by eq. (||), and its detailed form need not be precised. Without background, the free energy T' for the set of 
particles with pair potential Vp can be formally expanded in a Mayer diagrammatic representation , with a leading 
term (1/2) / n{r)vp{r, Y')n{r')cfr (Pr' for the excess (over ideal) part of J^'. The presence of a neutralizing background 
cancels this mean-field electrostatic term and the intrinsic free energy functional of the inhomogeneous OCP becomes: 



T[n{r)] = T'[n{r)] n{r) Vp{r,r') n{r') (fr d 

The local chemical potential reads 



2r'. 



Mr) = 
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and the second functional derivative of J- yields: 
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where the variations of the excess contribution to give rise to the usual direct correlation function 
short-range part given by 



(7) 

(8) 

(9) 
(10) 

(11) 
having a 



c,,(r,r') =c(r,r') + pVpir,r'). (12) 
Note that the chemical potential of the OCP on the sphere coincides with /i(0) for the optimum density profile (||). 
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Equation ( p| ) emphasizes the short range dependence of the chemical potential on a density perturbation. Conse- 
quently, fi{0) is the same in a finite TV-particle OCP in the central potential (r) and in the limit iV — » oo with an 
external potential V^{r) ensuring the same density variation around the origin as expression (||), namely: 

n{r) = n,(^l~^^ + ... (13) 

For the purpose of the present analysis, it is sufficient to truncate (H) after second order in r, as becomes clear below. 
The knowledge of the finite-size correction to the chemical potential for the OCP on the sphere then amounts to 
computing the shift Sfi{0) induced by switching V^{r) starting from the infinite homogeneous planar OCP with 
density (corresponding to the stereographic projection of the "spherical" plasma in the thermodynamic limit 
R — > oo). The density variation caused by the addition of Vj^{r) reads Sn{r) ~ —nsr'^/{2E?) and induces the shift 



5n{v)d?v (14) 



n(r) 

c.JOr^d^r, (15) 



where the direct correlation function to be considered is that of the homogeneous reference planar OCP. From the 
sum rule (^), we finally obtain: 



and eq. (0) is recovered. 
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FIG. 1. Stereographic projection from the North pole onto the plane V. A point Af on the sphere is projected onto P, with 
Cartesian coordinates r (r = at the South pole 5"). 
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